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ABSTRACT 

We present adaptive optics-assisted integral field spectroscopy around the Ha or H/3 
lines of 12 gravitationally lensed galaxies obtained with VLT/SINFONI, Keck/OSIRIS 
and Gemini/NIFS. We combine these data with previous observations and investigate 
the dynamics and star formation properties of 17 lensed galaxies at 1 < z < 4. 
Thanks to gravitational magnification of 1.4 — 90 x by foreground clusters, effective 
spatial resolutions of 40—700 pc are achieved. The magnification also allows us to probe 
lower star formation rates and stellar masses than unlensed samples; our target galaxies 
feature dust-corrected SFRs derived from Ha or H/3 emission of ~ 0.8 — 40M©yr _1 , 
and stellar masses M* ~ 4 x 10 8 — 6 x 10 10 M©. All of the galaxies have velocity 
gradients, with 59% consistent with being rotating discs and a likely merger fraction 
of 29%, with the remaining 12% classed as ‘undetermined.’ We extract 50 star-forming 
clumps with sizes in the range 60 pc - 1 kpc from the Ha (or H/3 ) maps, and find that 
their surface brightnesses, H c i ump and their characteristic luminosities, Lq, evolve to 
higher luminosities with redshift. We show that this evolution can be described by 
fragmentation on larger scales in gas-rich discs, and is likely to be driven by evolving 
gas fractions. 
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1 INTRODUCTION 


Advances in observational facilities and instrumentation 
over recent decades have led to a rapid accumulation of data 
on the statistical properties of galaxy populations and their 
evolution with cosmic time. Most notably, it is established 
that the cosmic star formation rate density peaked around 
1 and has sin c e declined by an order of magnitude 


(iMadau et al.l 19 961: Lilly et al.1 Il996l : iHopkins &; Beacoml 
20061 : ISobral et al.1 l2012h , and that more than half of the 
present-day stellar mass had been for med by z ~ 1 (e.g. 
iDickinson et al.l [20031 : IPatel et ai1l2013h . 
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The challenge now is to interpret these statistical trends 
and to uncover the underlying physical processes at work. 
How did these early galaxies form their stars, and what con¬ 
ditions have changed over time to result in the observed 
evolution of their integrated properties? 

Star formation is a complex process driven by the inter¬ 
play of competing effects: self-gravitating gas is supported 
against collapse by turbulent motion and by feedback in the 
form of radiation pressure from hot, young stars and active 
galactic nuclei. An understanding of the evolution of star¬ 
forming galaxies therefore requires an understanding of the 
evolving dynamics of galaxies at different epochs. A popula¬ 
tion of rapidly star-forming galaxies at high redshift implies 
the presence of large gas reservoirs. Since the dynamics of 
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gas and stars are intrinsically different - stars comprise a 
collisionless system, while gas is collisional and can dissi¬ 
pate energy - it is to be expected that high-redshift galaxies 
will have systematically different dynamics than those found 
locally. 

The advent of near-infrared tntegral Field Units (IFUs) 
on 8-10m class telescopes presented the opportunity to study 
two-dimensional velocity fields. At z > 1, the optical nebu¬ 
lar emission lines such as Ha are redshifted into the near- 
infrared, with [On] visible out to z ~ 5, providing bright 
tracers of the underlying dynamics in ionised gas across a 
large redshift range. 

IFU studies of the dynamics of high -z galaxies quickly 
revealed that a significant fraction, around 1/3, are rotat¬ 
ing systems in place as early as z ~ 2, but that they 
show higher velocity dispersions than local galaxies. They 
also found that ~ 1/3 are ‘dispersion dominated’ non¬ 


systems i Forster Schreiber et al. 2006b; 

Genzel et al. 2006: 

Law et al. 2007|, 2009; Wright et al. 20071. 

2009: Eoinat et al.l 

2009. 2012. see iGlazebrook 2013 for a 

comprehensive re- 


view) . 

A common feature in Ha maps of high-redshift galax¬ 
ies from IFU studies is the irregular, clumpy morphol¬ 
ogy of the Ha emission. These have also been observed in 
high-resolution imaging (Cowie et al. 1995; Elmegreen et alJ 


120041 : lEhnegreen fc Elmegre en 2005; Elmegreen et alJ 20091 ) 

and termed ‘clump cluster’ or ‘chain’ galaxies. The prevail¬ 
ing view is that these clumps form from internal gravita¬ 
tional instabilities in gas-rich discs (lElmegreen~et aDM, 


l200d : lGenzel et, aflteOOSltBournaud et, al.l201(t l. Their ubiq- 
uity in galaxies with ordered rotation supports this view, as 
clumps forming through major mergers would be ex pected 
to disrupt the dynam i cs (e.g. Bournaud &; Elmegreenl 12009 : 
iBournaud et al.ll201ll : iDekel et al.ll20od b 

The study of dynamics and star formation morpholo¬ 
gies of high-redshift galaxies is hampered by spatial resolu¬ 
tion. Even with adaptive optics, the current generation of 
telescopes can achieve resolution of ~ 1 — 1.5 kpc at z ~ 2. 
This limits the number of spatial resolution elements, as not 
only are galaxies intrinsically smaller at high redshift, but 
cosmological surface brightness dimming oc (1 + z) 4 limits 
observations to the bright central regions. Beam-smearing 
also affects measurements of velocity dispersion, which in¬ 
cludes contributions from the underlying velocity gradient 
as well as local turbulent motion, and can smooth out ro¬ 
tation curves to give slowly rotating g alaxies the appear¬ 
ance of being dispers ion-dominated (e.g. IWright et al] 
iNewman et alJl2013h . 


1.1 Gravitational magnification 

Enhanced spatial resolution can be achieved with current 
facilities by targeting galaxies that benefit from strong grav¬ 
itational lensing by massive foreground clusters. As well as 
stretching the galaxy images and thus increasing the spatial 
resolution of observations, lensing conserves surface bright¬ 
ness, which means that the total flux is effectively magnified. 
This enables us to study galaxies whose intrinsic luminosi¬ 
ties lie below the detection limits of current surveys, opening 
up more ‘normal’ galaxies, where unlensed surveys have by 


necessity focused on the more extreme star-forming popula¬ 
tion. 

By taking advantage of magnification due to strong 
gravitational lensing, high-resolution dynamics in high- 
redshift galaxies have been obs erved, revealing a high frac¬ 
tion to have ordered r otation (INesvadba et al.l 20061. 2007; 


Swinbank et al.l 20061: IStark et al.l 120081 : I Swinbank et al l 
2009: 1 Jones et al l2010h . Star-forming clumps are also visible 


in these galaxies, and due to the high spatial resolution can 
be observed on ~ 100 pc scales, enabling direct comparisons 
with the Hll regions in which stars form in t he local Universe 
( Swinbank et al.l 120091 : I Jones et"ahl l2010h . I Livermore et al.l 
( 2012ah used Ha narrowband imaging to study clumps in 
star-forming galaxies at 0 < z < 2 and showed that both 
their surface brightnesses and luminosity functions evolve 
with redshift, with the higher-redshift galaxies having more 
massive, brighter clumps with higher surface brightnesses. 
The high-z clumps can be explained by the same forma¬ 
tion process as Hll regions in local galaxies: a marginally 
stable disc fragments on scales related to the disc’s Jeans 
mass. At high redshift, higher gas fractions cause collapse 
on larger scales, leading to clumps large enough to dominate 
the galaxy’s morphology. This factor alone, however, over ¬ 
predict s the luminosities of high- z clumps. I Livermore et al.l 
(2012al ) therefore suggested a contribution from the epicyclic 
frequency, which is higher at high redshift and acts to sta¬ 
bilise the disc, causing collapse on smaller scales. 

This model effectively explains the evolution in ob¬ 
served clump properties, but relies on the dynamics of the 
galaxies, which were derived from disc scaling relations and 
not measured in the data. In this paper, therefore, we use 
integral held spectroscopy of lensed galaxies to examine the 
evolution of galaxy dynamics in combination with the prop¬ 
erties of star-forming clumps. The paper is organised as fol¬ 
lows: we present the sample and describe the data reduction 
and derivation of galaxy properties in Section [2] In Section 
02 we analyse the results and discuss first the dynamics and 
then the star-forming clumps, and we discuss the formation 
and evolution of star-forming clumps in Section [H Finally, 
we summarise our conclusions in Section [5j Throughout, 
we adopt a ACDM cosmology with Ho = 70kms -1 Mpc -1 , 
Ua — 0.7 and U m — 0.3. 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Integral Field Spectroscopy 

Table [T| lists the galaxies in our sample, which were se¬ 
lected from clusters with existing mass models, and to have 
spectroscopically confirmed redshifts that place their Ha or 
H /3 emission in windows of atmospheric transmission in the 
near-infrared and away from OH emission lines. In order to 
ensure detection in ~ 5 hours, we targeted galaxies with in¬ 
tegrated emission line fluxes of > 10“ 16 ergs -1 cm -2 , known 
either from prior s pectrosc opy or from spectroscopic pre¬ 
imaging (e.g. iRichard et alJl201lh . 

The new sample comprises twelve galaxies, of which 
ten were observed with t he SINFONI Integral F ield Unit 
(IFU) on the ESO/VLT (|Eisenhauer et all 120031 ) between 
2009 April 30 and 2012 July 16. The SINFONI targets were 
selected to be close to sufficiently bright stars for the use 
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Table 1. Targets and observations. 


Name 

(a) 

02000 

0 b) 

<52000 

(b) 

z 

0>) 

Instrument 

Filter 

Emission line 

Exposure time (ks) 

MACS0744-system3 

07:44:50.90 

+39:27:34.5 

1.28 

NIFS 

H 

Ha 

19.2 

MACS1149-arcAl.l 

11:49:35.30 

+22:23:45.8 

1.49 

OSIRIS 

Hn3 

Ha 

17.1 

MACS0451-system7 

04:51:57.22 

+00:06:21.3 

2.01 

SINFONI 

H 

H/3 

19.2 

A1413-arc2.1a 

11:55:18.31 

+23:23:55.5 

2.04 

SINFONI 

K 

Ha 

21.6 

A1413-arc2.1b 

11:55:18.31 

+23:23:55.5 

2.04 

SINFONI 

K 

Ha 

21.6 

A1835-arc7.1 

14:01:00.99 

+02:52:23.2 

2.07 

SINFONI 

K 

Ha 

19.2 

MS1621+26-systeml 

16:23:35.56 

+26:34:29.1 

2.14 

SINFONI 

K 

Ha 

21.6 

RXJ1720+26-arcl.l+1.2 

17:20:10.26 

+26:37:27.4 

2.22 

SINFONI 

K 

Ha 

24.0 

A1689-arc2.1 

13:11:26.52 

-01:19:55.5 

2.54 

SINFONI 

H+K 

H/l 

19.2 

A1689-arcl.2 

13:11:26.43 

-01:19:57.0 

3.04 

SINFONI 

H+K 

H/3 

19.2 

A2895-arcl. 1+1.2 

01:18:11.15 

-26:58:03.9 

3.40 

SINFONI 

K 

H/3 

19.2 

A2895-arc2.2 

01:18:10.45 

-26:58:14.2 

3.72 

SINFONI 

K 

H/l 

26.4 

C10024-arcl.l 

00:26:34.42 

+17:09:55.4 

1.68 

OSIRIS 

Hn5 

Ha 

16.5 

C10949-arcl 

09:52:49.78 

+51:52:43.7 

2.39 

OSIRIS 

Hn4 

Ha 

19.2 

MACS0712-systeml 

07:12:17.51 

+59:32:16.3 

2.65 

OSIRIS 

Kc5 

H/3 

16.2 

MACSJ0744-arcl 

07:44:47.82 

+39:27:25.7 

2.21 

OSIRIS 

Kn2 

Ha 

14.4 

CosmicEye 

21:35:12.73 

-01:01:43.0 

3.07 

OSIRIS 

Knl 

H/3 

21.6 


Notes: (a) To identify the arcs, we use the cluster names and the arc ID from the published lens model or the system number where 
we observe multiple images of the same galaxy, except for the Cosmic Eye which is also known as LBG J213512.73-010143. MACS 
clusters are truncations of MACSJ0744.8+3927, MACSJ 1149.5+2223, MACSJ0451.9+0006, MACSJ0712.3+5931 and MACSJ2135.20102 
(lEbeling et al.lliooll , 120071 , l2010l : iMann fc EbelineJl2012lV (6) Positions and redshifts are given for the target lensed arcs. 


of NGS+AO, resulting in a median FWHM = 0.2”. As the 
lensed arcs are extended, usually over several arcseconds, 
we used the 8” x8” field of view with a spaxel size of 0.25”. 
Due to the elongated shape typical of lensed arcs, the tar¬ 
gets were kept in one half of the IFU and nodded across in 
an ABBA sequence. We observed each target for six ABBA 
sequences in the H or iC-band filter according to the red- 
shifted position of the target emission line (see Table [T}. We 
also observe [Nil] in the Ha targets and [Om] in the H [3 
targets; the properties of these lines will be discussed in a 
future paper. 

The data were reduced using the esorex package, 
which performs flat-fielding and wavelength calibration and 
reforms the image into a data cube with two spatial and 
one spectral dimension. It also carries out sky subtraction 
by subtracting each B frame from its closest A frame, which 
with our observing strategy results in cubes that contain 
two images of the target, one positive and one negative. 
Standard stars were observed on the same nights and in the 
same filters as the science exposures, and were reduced in 
the same way. We used the standard stars to individually 
flux-calibrate each cube before combining them. 

To combine the cubes, we first cut them in half to sepa¬ 
rate the two images of the target, and subtracted the nega¬ 
tive image from the positive one. We then aligned the cubes 
by collapsing them into continuum images. As the targets 
all lie behind cluster lenses, they commonly have elliptical 
cluster members close to their line of sight. Where possible, 
a foreground elliptical galaxy was intentionally positioned 
within the SINFONI field of view to allow careful alignment 
between cubes. Where this was not possible (in the case of 
the two arcs lensed by the cluster Abell 2895), we aligned 
the cubes by centring on a bright feature within the lensed 
galaxy. Once aligned, the cubes were median-combined with 


the IRAF task imcombine with the crreject algorithm for 
cosmic ray rejection. 

In addition to the ten targets observed with SIN¬ 
FONI, we observed a z = 1.28 galaxy lensed by the clus¬ 
ter MACS J0744+3927 with Gemini/ Near-Infrared Inte¬ 
gral Field Spectrometer (NIFS: [McGregor et al.ll2003h . The 
lensed target is denoted MACS0744-system3. As the target 
is too long to fit in the 3” x 3” field of view (see first panel 
of Figure [I]), we observed in an ABC sequence, where two 
halves of the arc were positioned diagonally across the field 
in the A and C frames, and the B frame was a blank field 
used for sky subtraction. The A and C frames were chosen 
to overlap at the foreground elliptical in the middle of the 
arc to enable precise alignment of the A and C frames. We 
observed for 8 ABC sequences with the 77-band filter. Data 
reduction was carried out using the gemini package in IRAF 
of both the science frames and standard stars observed on 
the same nights and with the same setup as the science ob¬ 
servations. Each cube was then individually flux-calibrated 
and combined as for the SINFONI observations described 
above. 

One further target, the spiral galaxy at z = 1.49 lensed 
by the cluster MACS J1149.5+2223 (MACS1149-arcAl.l), 
was observed with Keck/OSIRIR T he o b serva tions and 
data reduction are described bv lYuan et al.1 ( 2011 5. who use 
the data to show that the galaxy’s metallicity gradient is 
strongly negative. 

Throughout this paper we also make use of the sam¬ 
ple of six le n sed galaxies observed with Keck/OSIRIS from 
Ijones et al.l (2010). They presented the resolved kinemat¬ 
ics of these galaxies, which we adopt in the first half of 
this paper for comparison with our new data. They also 
demonstrated that individual star-forming clumps could be 
extracted from the data and discussed their properties in re- 
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lation to local Hu regions, deriving luminosity densities up 
to 100 x higher than local star-forming regions. So that they 
can be combined with the new data in a self-consistent man¬ 
ner, we undertake a new analysis of the clump properties in 
this sample in the latter half of this paper, but ou r results 
are en tirely consistent with those presented bv ljones et al.l 
(|2010h . One of these galax i es - M ACS0451-system7, ob¬ 
served in Ha bv ljones et al.l ( 2010 ) - is also included in the 
SINFONI sample, where we observed it in H (3 and [Om]. The 
combined sample of lensed arcs thus comprises 17 galaxies 
at 1.28 < * < 3.72. 


2.1.1 Unlensed samples 

For comparison purposes, we refer throughout this paper to 
three other IFU surveys of high-redshift galaxies, which are 
not lensed. For ease of reference, we summarise the main 
properties of these surveys as follows: 


SINS+AO: This is a subset of the Spectroscopic Imag¬ 
ing survey in the Ne ar-infra red with SINFONI (SINS; 
iForster Schreiber et al.l l2006a l) sample whi ch ha s been ob¬ 
served with adaptive optics " ^Newman et al.ll2013l ). The sam¬ 
ple comprises 38 star-forming galaxies selected via their 
UV or optical emission, and has z ~ 2.2. These galax¬ 
ies sample the massive end of the star-forming ‘main se¬ 
quence,’ with stellar masses of 10 9 ' 2 — 10 11 ' 5 M© and SFRs 
of 13 — 850M©yr -1 . The spatial resolution achieved with 
SINFONI+AO is ~ O'/2, equivalent to 1.7 kpc. 

SHiZELS: A sample of nine galaxies at 0.84 < z < 2.23 
drawn from the High-Z Emission Line Su rve y (HiZELS), 
observed with SINFONI+AO (|Swinbank et al.l 1 20121 ). The 
spatial resolution of O'/I equates to 0.8kpc at the sample’s 
median redshift, z = 1.47. The sample covers a similar range 
of stellar masses (M* ~ 10 9 — 10 11 M©) to the SINS sample, 
but with lower SFRs of 7 d= 2M©yr -1 . 

WiggleZ: This sample comprises 13 z ~ 1.3 galaxies se¬ 
lected from the WiggleZ Dar k Energy Survey and observe d 
with Keck/OSIRIS+LGSAO ([Wisnioski et al.ll201lLl2012ah . 


The spatial resolution of these observations is ~ O'/I, or 
0.8 kpc. The range of stellar masses covered is similar to the 
other two surveys at 10 9 ' 8 M© < M* < 10 11 ' 7 M©, and SFRs 
are in the range 30 — 200M©yr -1 . 


2.2 HST imaging and lens modelling 

HST colour images (where available) of the target galaxies 
are shown in the top-left panels of Figure [T] In order to de¬ 
termine the intrinsic properties of the target galaxies, it is 
necessary to account for the effects of lensing. To do this, 
we use the best-fit cluster mass models, the details of which 
are given in Table [2] Where the cluster’s critical line at the 
redshift of the lensed arc crosses the field, it is overlaid in 
red. The critical line denotes regions of theoretically infinite 
magnification, and arcs which cross this line are multiply im¬ 
aged within the IFU field of view. To reconstruct the source- 
plane images of t he galaxies, we use lenstool ([Knei 
I Julio et aLll2007l ) to ray-trace each pixel to its source-plane 
origin. After using a central pixel to obtain the position of 
the galaxy in the source plane, we construct a regular grid 
in the source plane, ray-trace each pixel to the image plane 
and obtain the spectrum in that position by interpolating 


between pixels in the IFU data cube. In the cases of A1835- 
arc7.1 and A2895-arc2.2, the high magnification gradient in 
the vicinity of the critical line causes this method to omit 
data from the image plane. For these galaxies, we therefore 
carry out the process in reverse by ray-tracing each pixel 
from the IFU data cube to the source plane. This results in 
an irregular grid in the source plane, which is gridded into 
square pixels using a delauney tesselation. Once the source 
plane cubes have been constructed, we apply conservation 
of surface brightness to each pixel to obtain the instrinsic 
source plane flux. The magnification factor [i for the galaxy 
given in Table [I] is then the ratio of image- to source-plane 
flux measured from the H a or H/3 emission lines. 

With both methods, the source plane pixel scale is cho¬ 
sen so that each pixel in the IFU data cube is represented 
in the source plane, and we impose a lower limit of 0.001” 
to ensure manageable file sizes. The pixel scale is therefore 
dictated by the direction or region of highest magnification. 
Gravitational lensing usually acts preferentially in one direc¬ 
tion, and galaxies lying close to the critical line will experi¬ 
ence strong magnification gradients. Hence, the source plane 
cubes are oversampled both in the direction of lower mag¬ 
nification and in regions of the galaxy that lie furthest from 
the critical line. To estimate the actual resolution achieved, 
we apply the same reconstruction as described above to the 
standard star observations and measure the FWHM in the 
image and source plane. The result is an ellipse, shown in 
magenta in Figure [I] We remind the reader that while lens¬ 
ing brings a major improvement in the spatial resolution 
that can be achieved, this gain is not without its problems. 
For example, in many systems the magnification is much 
higher in one direction than the other, leading to a distorted 
image of the galaxy. While we find that reconstruction of the 
original image is not problematic for the lenses considered 
here, this limitation should be born in mind when assessing 
the overall morphology of the source. The effective PSFs 
shown in Figure [I] illustrate the extent of anisotropy in each 
system, and we give values for the magnification in each 
direction in Table [2] 

In order to estimate the errors on the magnification 
factors for each galaxy, we use the family of 100 best-fit 
lens models. We use lenstool to reconstruct the Ha or H/3 
intensity maps for each galaxy for each one of the possi¬ 
ble lens models, and then measure the ratio of image- to 
source-plane flux in each. The la deviation in the derived 
magnifications are given in Table |T| as the error on the to¬ 
tal magnification. As we measure the magnification for the 
purpose of this paper from the Ha (or H/3) line intensity 
(effectively a weighted mean), the largest errors are found 
where there is a large magnification gradient across the arc 
without strong constraints on the mass distribution of the 
primary lensing source. The largest fractional error, 35%, is 
found in MACS0744-system3, due to the strong magnifica¬ 
tion gradient across the image arising from lensing by the 
foreground cluster galaxy positioned along the line of sight 
to the arc. The smallest error, 5%, is found in A1835-arc7.1, 
which also has a strong magnification gradient, but in this 
case the two images straddling the critical line provide good 
constraints on the precise position of the line. We note that 
changing the lens model tends to result in the galaxy image 
being stretched differently, so that the effective resolution 
described above would change, but the underlying morphol- 
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Figure 1. Image plane colour images (where available) and reconstructed source plane maps of the target galaxies. The top-left image 
in each panel shows an HST image of the arc constructed from archival ACS and/or WFC3 data, in colour where there are multiple 
bands in the archive. The filters used are given in the lower-left corner of the image. The critical line at the target galaxy’s redshift is 
overlaid in red. Where the IFU field of view cannot cover the entire arc, it is shown as a dashed white box. The lower-left image shows 
the Ha or H (3 emission line flux, aligned to the same astrometry as the HST image for context. The right four images in each panel 
show the source plane reconstructions: from top left to bottom right, they show the Ha or H j3 intensity maps, velocity fields, line-of-sight 
velocity dispersion and velocity field residuals after subtracting the best-fit disc model. The star-forming clumps are contoured over the 
intensity maps. Contours overlaid on the velocity field show the best-fit disc model after smoothing by the effective source plane PSF, 
shown by a magenta ellipse. As surface brightness is conserved by lensing, the image and source plane intensity maps are displayed in 
terms of surface brightness, with the same scaling. 


ogy cannot be materially altered. Crucially for the discussion 
in the latter half of this paper, new features such as clumps 
cannot be created during the lensing recontruction. 


2.3 Dynamical maps 

With both the image- and source-plane cubes, we fit the in¬ 
tensity, velocity and velocity dispersion of the emission lines 
in each pixel using a y 2 minimisation technique. In each 
pixel, we simultaneously fit the Ha and [Nn]A6583 or H/3 
and [Oiii]AA4959,5007 lines. To reduce the number of pa¬ 
rameters in the fit, we require all lines to be of the same 
velocity and velocity dis persion, and we impose a ratio of 
[Oiii]A5007/A 4959 = 3 (|Storev k, Zeipper3l2000h . The fit is 
accepted if it results in a Ay 2 > 25 compared to a straight- 
line fit, equivalent to a signal-to-noise of S/N > 5. If no 
fit is obtained, we adaptively bin up to 3x3 pixels in the 
image plane, or the equivalent area of 3x3 PSF areas in 
the source plane. We deconvolve the velocity dispersion for 
spectral resolution by subtracting in quadrature the median 


width of sky lines measured from blank exposures taken with 
the standard stars. The high spatial resolution of the source- 
plane cube means that most of the contribution to the ve¬ 
locity dispersion of the velocity gradient across the galaxy is 
removed. To remove the remainder, we subtract in quadra¬ 
ture the velocity gradient across each pixel, measured over 
the PSF. For each parameter, we also estimate the errors 
in each pixel. Taking the best ht to the emission lines, we 
vary each parameter in turn, while allowing the others to 
find their new minima, until we obtain a Ay 2 ^ 1. 

We show the resulting maps in Figure [I] In the bottom 
left of each panel is the Ha or H (3 emission line intensity 
from the image-plane IFU data, mapped to the same as¬ 
trometry as the HST image. The source plane reconstruc¬ 
tions are in the right-four images of each panel, showing the 
Ha or H/3 emission line intensity, the velocity and velocity 
dispersion maps and the velocity residuals after subtracting 
the best-fit disc models (contoured over the velocity map 
and described in Section ETT1) . We note that as the velocity 
and velocity dispersion are fixed between multiple emission 
lines, their values are dominated by the line with the high- 
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Figure \l\ (continued) 
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Figure \l\ (continued) 
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Figure \l\ (continued) 


Table 2. Gravitational lens properties of the sample. 


Name 

V 

Win 

Resolution (pc) 

Lens model reference 


(a) 

C b ) 

(c) 


MACS0744-system3 

63 d= 22 

21 x 3 

134 ± 47 

Jones et al. (2010): Korneut et al. (2011) 

MACS 1149-arcAl. 1 

16 ± 2 

5x3 

169 ± 21 

Smith et al. (2009) 

MACS0451-system7 

90 ± 7 

63 x 1 

84 ± 7 

Jones et al. (2010) 

A1413-arc2.1a 

24 ± 3 

18 x 1 

299 d= 37 

Richard et al. (2010) 

A1413-arc2.1b 

29 ± 2 

21 x 1 

265 d= 18 

Richard et al. (2010) 

A1835-arc7.1 

63 ± 3 

37 x 2 

118 d= 6 

Richard et al. (2010) 

MS1621+26-systeml 

8 ± 2 

4x2 

722 ± 165 

Richard et al. (in prep) 

RXJ1720+26-arcl.l+1.2 

86 d= 18 

14 x 6 

560 d= 115 

Richard et al. (2010) 

A1689-arc2.1 

15 ± 3 

21 x 1 

96 d= 19 

Limousin et al. (2007) 

A1689-arcl.2 

61 ± 11 

25 x 2 

249 =b 47 

Limousin et al. (2007) 

A2895-arcl. 1+1.2 

7 ± 1 

10 x 1 

454 ± 88 

May et al. (in prep) 

A2895-arc2.2 

9 ± 2 

7 x 1 

167 d= 40 

May et al. (in prep) 

C10024-arcl.l 

1 d= 0 

2x1 

508 ± 55 

Jones et al. (2010) 

C10949-arcl 

7 ± 2 

6x1 

138 ± 38 

Jones et al. (2010) 

MACS0712-systeml 

28 ± 8 

19 x 1 

40 ± 11 

Jones et al. (2010) 

MACSJ0744-arcl 

16 d= 3 

9x2 

83 ± 15 

Jones et al. (2010) 

CosmicEye 

28 ± 3 

8x3 

96 ± 10 

Jones et al. (2010) 


Notes: (a) fi is the total magnification of the source, derived from the ratio of the emission line flux in the image plane to that in the 
source plane. fiu n gives the linear magnification factors in the direction of greatest magnification and perpendicular to that axis, (c) 
Resolution is given in the direction of greatest magnification. 
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Table 3. Dynamical properties of the sample 


Name 

Inclination 

e 

PA 

V 2.2 
(km s -1 ) 

a 

(km s _1 ) 

v/r 

(km s _1 kpc _1 ) 

v/a 

Disc/ Merger/ 
Undetermined 

MACS0744-system3 

61 

266 

80 d= 10 

60 ± 20 

230 d= 30 

1.4 d= 0.8 

Disc 

MACS 1149-arcAl. 1 

45 

86 

59 ± 3 

50 ± 10 

88 ± 4 

2 ± 1 

Disc 

MACS0451-system7 

55 

100 

100 ± 10 

77 ± 9 

220 ± 30 

1.6 ± 0.6 

Disc 

A1413-arc2.1a 

50 

172 

17 ± 4 

50 ± 10 

24 ± 5 

0.5 d= 0.3 

Disc 

A1413-arc2.1b 

66 

249 

20 ± 3 

60 ± 20 

8 ± 1 

0.6 ± 0.3 

Merger 

A1835-arc7.1 

89 

236 

40 ± 10 

70 ± 10 

4 ± 1 

0.2 ± 0.2 

Undetermined 

MS1621-systeml 

84 

247 

126 ± 4 

60 ± 20 

47 ± 1 

3 d= 1 

Merger 

RXJ1720-arcl.l+1.2 

70 

144 

63.2 ± 0.5 

60 =b 20 

27.1 ± 0.2 

0.4 ± 0.2 

Merger 

A1689-arc2.1 

81 

38 

64 ± 7 

55 ± 1 

180 ± 20 

1 ± 1 

Disc 

A1689-arcl.2 

80 

234 

78 ± 4 

90 ± 10 

53 ± 4 

1.3 ± 0.4 

Merger 

A2895-arcl.l+1.2 

85 

260 

60 ± 10 

60 db 20 

62 ± 3 

1.9 d= 0.7 

Undetermined 

A2895-arc2.2 

47 

276 

200 ± 10 

60 ± 20 

530 ± 30 

4 ± 2 

Disc 

C10024-arcl.l 

50 

122 

91 d= 8 

63 ± 6 

15 ± 3 

2.6 d= 0.3 

Disc 

C10949-arcl 




68 d= 7 

16 ± 3 


Merger 

MACS0712-systeml 

40 

351 

33 d= 3 

75 ± 7 

25 ± 3 

0.81 ± 0.08 

Disc 

MACSJ0744-arcl 

45 

124 

118 d= 9 

110 d= 10 

210 d= 20 

2.4 d= 0.2 

Disc 

CosmicEye 

55 

113 

51 d= 4 

44 d= 4 

55 ± 5 

1.9 ± 0.2 

Disc 


Notes: Inclination and PA are derived from the best-fit disc models. V 2.2 is the velocity at 2.2 r 1 / 2 , where r 1 / 2 is the half-light radius. 
Velocity dispersion, cr, is the luminosity-weighted mean l ocal value, after de convolving for local velocity gradient. The velocity gradient 
v/r is measured at r 1 / 2 . Targets below the line are from I Jones et al ] 1 I 20 KJ) with values reported in their Table 2, except for V 2 . 2 , v/r 
and v/a, which we measure for consistency with the new data. No values for 6 , PA or V 2.2 are given for C10949-arcl as this galaxy cannot 
be fit by a disc model. 


est signal-to-noise. In the case of the H/3/ [Om] observations, 
the [Om]A5007 line is far brighter than the H f3 line; in some 
cases we thus obtain velocity and velocity dispersion values 
in pixels where no H/3 is measured. We n ote that the method 
above is the same as that used bv I Jones et al.l ( 2010 ), so the 
properties derived from the two samples are self-consistent. 


2.4 Disc modelling 

To each of the velocity fields, we follow Ijones et~al . (l2010h 
and attempt to fit a rotating disc model and hence infer 
the true disc rotation speed and inclination. The model is 
described by six parameters: the disc centre x and y, the 
asymptotic rotational velocity I’asym, the turnover radius r*t, 
the position angle </> and inclination 0. The values of these 
parameters are constrained so that x and y are within the 
range of the data, and rt < 5 0kpc. Rotation is described by 
an arctan function (|Courteaulll997h of the form 


v(r) — 



rasym arctan 



(i) 


The fit is carried out with an iterative procedure 
whereby 10, 000 random sets of parameters are generated, 
and the resulting 2D velocity models are smoothed by the 
effective source plane PSF of the galaxy in order to account 
for the effect of asymmetric magnification. The % 2 is then 
calculated for each one. The parameter space is contracted 
by discarding up to 10% of the worst fitting models in each 
iteration. This process is repeated up to 100 times, and the 
fit is deemed to have converged once all models in a gen¬ 
eration have A% 2 < 1. The smoothed best-fit disc models 
are contoured on the observed velocity fields in Figure [l] 
and the velocity residuals after subtracting the best-fit disc 
model are shown in the bottom-right image of each panel. 


The dynamical properties of the sample are given in Table 

El 

While it is always possible to fit a rotating disc model 
to the data, this is insufficient to determine that the galaxy 
is a rotating disc. A rotating disc should have symmetric ve¬ 
locity and velocity dispersion fields, the latter peaking in the 
centre of the galaxy. The traditional requirement that the 
velocity field form a ‘spider’ pattern is rarely found in lensed 
galaxies due to the asymmetric magnification (the contours 
overlaid in Figure [T] demonstrate how much smoothing by 
the elliptical source plane PSFs distorts even an idealised 
rotation field ). 

Ijones et al. liol3) illustrate the symmetry of the kine¬ 
matics in their sample by extracting one-dimensional profiles 
along the kinematic axis, concluding that 5/6 of the galaxies 
have velocity profiles indicative of rotation, with one likely 
merger (C10949-arcl). In order to examine the symmetry 
in our new data, we extract a one-dimensional profile of 
the velocity and velocity dispersion from each galaxy. We 
do this by using the best-fit disc model to identify the dy¬ 
namical centre and major kinematic axis of the galaxy. The 
velocity and velocity dispersion are extracted along a slit 
five pixels wide, and binned into independent resolution ele¬ 
ments along the long axis according to the galaxies’ effective 
source plane PSF. The resulting dynamical profiles for the 
12 galaxies in our sample are shown in Figure [2] The best 
fit one-dimensional rotation curves and exponential velocity 
dispersion profiles (where applicable) are overlaid in red. 

All of the galaxies in the sample exhibit velocity gra¬ 
dients, and all can be fit by an arctan function. However, 
we only see indications of a turnover in the rotation curve 
in 9/12 of the sample. Furthermore, the velocity dispersion 
profiles are less ordered. In Figure [2] we attempt to fit an 
exponential profile to the velocity dispersions. In 50% of 
cases the exponential function is a good fit to the data. In 
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Figure 2. One-dimensional rotation curves (left) and velocity dispersion profiles (right), extracted from a slit aligned along the axis of 
the best-fit disc model. For the galaxies whose dynamics are well fit by a rotating disc model, we overlay the best-fit rotation curve and 
velocity dispersion profile. 
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two cases, A1413-arc2.1b and RXJ1720-arcl.l+1.2, there is 
no central peak. A1413-arc2.1b is probably interacting with 
the larger A1413-arc2.1a, from which it is offset by only 
~ 1000kms -1 , and the morphology of RXJ1720-arcl.l+1.2 
with its extended tail suggests that it also may be an inter¬ 
acting system. In the remaining four galaxies - MACS0744- 
system3, A1413-arc2.1a, MS1621-systeml and A1689-arcl.2 
- there are suggestions of a double peak in the velocity dis¬ 
persion profile, which may also be indicative of a late-stage 
merger. Alternatively, the irregular and asymmetric profiles 
could be due to turbulence within the discs. 

The kinema tic prop e rties o f the galaxies are similar to 
those in the I Jones et al.l (2010) sample, which also exhib¬ 
ited rotation-like velocity fields with irregular velocity dis¬ 
persion profiles. The galaxies in t he n ew sam ple tend to 
be smaller than those from Ijones et al ] (|2010l b with me¬ 
dian r 1/2 ~ 0.5 kpc and 1.2 kpc respectively, but both lensed 
samples probe systematically smaller galaxies than unlensed 
studies such as SINS+AO, SHiZELS and WiggleZ, which 
have median r-i / 2 ^ 4.3 kpc, 2.3 kpc and 2.5 kpc respectively 

(INewman et aljfeoijl: ISwinbank et al.|[2012l: IWisnioski et alj 

l201lh . The two lensed samples are closely matched in ro¬ 
tational velocity with median rotational velocity measured 
at 2.2 r 1/2 of V 2.2 ~ 64 km s -1 and 51kms -1 for the new 
data and Ijones et~al . (l2010h sample respectively. The me¬ 
dian rotational velocity in the WiggleZ survey is marginally 
higher at 74kms _1 , with ~ 110 km s -1 in SHiZELS and 
~ 150kms _1 in the SINS AO sample. There are two fac¬ 
tors that bias lensing surveys towards more slowly rotating 
systems: firstly, that they are sensitiv e to smaller galaxies , 
which tend to be slower rotators (e.g. INewman et ai1l2013L 
see also Section [3.1.ID . and secondly, the high spatial reso¬ 
lution enables us to measure small velocity gradients that 
would be flattened by beam smearing in unlensed data and 
thus categorised as non-rotating. 

As the properties of the two lensed samples are similar, 
and differ systematically from the unlensed data, we discuss 
them in the remainder of this paper as a combined sample. 
In the interest of providing a comparison to other studies, 
we give our best estimate of the nature of the galaxies in 
Table [3J but caution that these are by no means certain. 


2.5 SED fitting and stellar mass estimates 

In order to estimate the extinction due to dust and the 
stellar masses of our sample, we model their spectral en¬ 
ergy distributions (SEDs) from available archival HST and 
Spitzer/IRAC imaging. We first subtract any foreground 
cluster members lying close to the line of sight to the tar¬ 
get galaxy, using the IRAF tasks ellipse and bmodel. We 
then degrade all of the imaging for each galaxy to the poor¬ 
est resolution, that of the longest-wavelength IRAC band. 
We then extract the photometry using an elliptical aperture 
large enough to encompass the galaxy, with a 2” annulus for 
sky subtraction. We apply aperture corrections to the IRAC 
fluxes based on the area of the aperture used. 

As the continuum light does not necessarily follow the 
same distribution as the nebular emission, we do not assume 
the same magnification factor from gravitational lensing; in¬ 
stead, the measured fluxes are corrected for lensing using the 
flux-weighted mean magnification within the aperture, from 
magnification maps created with LENSTOOL. The resulting 


magnification factors are within 1 a of the values given for 
the nebular emission in Table [T] The complete photometry, 
corrected for lensing, is given in Tables lAll IA2l and lA3l in the 
appendix. When calculating stellar masses, we omit the arcs 
in three clusters that do not have imaging in at least two 
bands either side of the 4000A X break (MS 1621, RXJ1720 
and Abell 2895), as well as the galaxy A1413-arc2.1b which 
is undetected in the HST imaging. 

We_perform the SED fitting using the CIGALE code 
(Noll et al. 1200 9) We use the stellar population models of 
iMarastonl ( 200 5). allow either exponentially decreasing or 
continuous star formation histories and constrain the old¬ 
est stellar populations to be younger than the age of the 
Universe at the target redshift. The CIGALE code creates 
SEDs from the far-ultraviolet to the infrared based on a 
dust-attenuated stellar population, infrared dust emission 
and spectral line templates. The best-fit galaxy parameters 
are derived with a Bayesian-like analysis from the distribu¬ 
tion of probability-weighted best-fit models. 

The stellar masses and dust extinction, Ay, obtained 
for each galaxy from the SED fit is given in Table 2] We 
find stellar masses of M* = 6 x 10 8 — 2 x 10 lo M©; thus, 
lensing allows us to probe systematically smaller galax¬ 
ies than un lensed samples su ch as SINS (median M* ~ 
3x 10 10 M^: lForster Schrei ber et al .ll2009h or SH iZELS (me¬ 
dian M* ~ 2 x 10 lo M©; ISwinbank et al.l [2012h . The dust 
extinctions we derive are in the range Ay = 0.4 — 1.1 with 
a median Ay ~ Q.9, similar to those of SINS (Ay ~ 0.8; 
Forster Schreiber et al. 12009 1 and SHiZELS (Ay ~ 0.9; 


Swinbank et al. 1 120121) . 


The stellar masses and dust extinctions of the 
Ijones et~al ] (I2010I) sampl e and A1835-a r c7.1 o f the new ob¬ 
servations are derived bv lRichard et all (|201lh and given in 
their Table 4. We note that the method used is similar to 
the one we employ, so the results can be directly compared. 


2.6 Integrated galaxy properties 

We now derive some integrated properties of the sample, so 
that we can relate the lensed galaxies to the population as 
a whole and explore how their internal structures evolve as 
a function of their global properties. 

A summary of the integrated properties of the sample 
is given in Table [ff] Total fluxes of the Ha or H /3 emission 
lines are calculated by summing every pixel in the source 
plane cubes with signal-to-noise > 5. We measure the total 
flux by fitting a Gaussian profile to the emission lines, and 
the ratio of the total flux in the image plane to that in the 
source plane gives the total magnification given in Table |T] 

Intrinsic star formation rates (SFRs) in Tableware cal¬ 
culated from the Hck (or H/3, assuming case B recombina¬ 
tion) flux, corrected for dust according to t he extinc tion Ay 
derived from SED fitting, by applying the lKennicu tt] (1998a) 
prescription. We correct this conversion to a Chabrier ini¬ 
tial mass function (IMF), which includes a turnover at the 
low-mass end and produces more realistic mass-to-light ra¬ 
tios. This correction has the effect of reducing the SFR by 
a factor 1.7x. 

To find the half-light radius ri/ 2 , we first construct an 
array in which each pixel value is the distance from the dy¬ 
namical centre, [z c , j c ], based on an estimated position angle, 
0, and inclination ( 0 ; the derivation of the dynamical centre, 











































Star formation rate (M 0 yr 


Resolved spectroscopy of gravitationally lensed galaxies 


15 


Table 4. Integrated properties of the sample 


Name 

Intrinsic fn a 

Intrinsic fn/3 

SFR 

U/2 

logM* 

Ay 


(10 _18 ergs 

-1 cm -2 ) 

(Meyr- 1 ) 

(kpc) 

(Mo) 


MACS0744-system3 

18 ± 2 


1.7 d= 0.2 

0.3 ± 0.1 

9.2 ± 0.4 

0.9 ± 0.4 

MAC S1149-arcAl.l 

12 ± 2 


1.2 ± 0.2 

0.6 ± 0.2 

9.3 ± 0.4 

0.5 ± 0.3 

MACS0451-system7 


5 ± 2 

5 ± 2 

0.18 ± 0.03 

8.8 ± 0.3 

0.7 ± 0.3 

A1413-arc2.1a 

34 ± 4 


6.6 ± 0.7 

0.46 ± 0.08 

8.9 ± 0.4 

0.4 ± 0.3 

A1413-arc2.1b 

10 ± 1 


3 ± 0.3 




A1835-arc7.1 

2.6 ± 0.3 


0.8 ± 0.2 

0.6 d= 0.2 

8.8 ± 0.2 a 

1.1 ± 0.3 a 

MS1621-systeml 

32 ± 4 


11 d= 1 

2 ± 1 



RXJ1720-arcl.l+1.2 

3.7 ± 0.5 


1.3 ± 0.2 

0.7 ± 0.3 



A1689-arc2.1 


0.7 ± 0.5 

1.1 ± 0.2 

0.3 ± 0.1 

9.5 ± 0.1 

0.9 :± 0.1 

A1689-arcl.2 


1 ± 0.5 

1.5 ± 0.3 

0.23 ± 0.08 

9.2 ± 0.2 

0.9 ± 0.2 

A2895-arcl.l+1.2 


10.4 d= 0.1 

22 ± 4 

0.5 ± 0.3 



A2895-arc2.2 


9 ± 3 

27 ± 5 

0.17 ± 0.09 



C10024-arcl.l 

130 d= 30 


27 ± 6 

1.8 d= 0.2 

10.4 ± 0.1 a 

1.1 ± 0.2 a 

C10949-arcl 

50 ± 10 


20 ± 6 

3.5 ± 0.9 

10.2 ± 0.5 a 

1.0 ± 0.0 a 

MACS0712-systeml 


2.8 ± 0.6 

5 ± 0.1 

0.5 ± 0.2 

10.5 ± 0.4 a 

0.8 ± 0.2 a 

MACSJ0744-arcl 

8 ± 2 


2.4 ± 0.5 

1.2 d= 0.2 

10.0 ± 0.1 a 

0.8 ± 0.2 a 

CosmicEye 


18 ± 2 

40 ± 5 

0.8 d= 0.2 

10.8 ± 0.1 a 

0.7 ± 0.1 a 


Notes: All quantities are corrected for lensing magnification. SFR is calculated from the H a or H/3 luminosity corrected for dust 
extinction based on A^/. The half-light radius r i/2 is based on the Ha or H/3 morphology as described in the text. A v and M* are 
est i mated f rom SED fitting of broadband photometry as described in the text, except for those marked (a), which are from 
[Ric hard et al. (2011). Targets below the line are from I Jones et al . (12010 s ). with SFRs given in their Table 2 and half-light radii reported 
bv lRichard et al.l (120111 s ). 



Figure 3. Total galaxy SFR against stellar mass for the lensed 
arcs for which we are able to derive stellar masses, in com¬ 
parison with other high-redshift kinematic studies. The lensed 
arcs are colour-coded by redshift. The dashed lines indicate 
the star-for ming galax y ‘main sequence’ at z = 0, 1 and 2.5 
(Whitaker [2C^|). Compari son data are from SHiZELS 

( SwjnbarfiTeriini2012h. WiggleZ JWisnioski et al l liolih and the 
SINS AO sample (iNewman et al.ll2013h . and the grey cross indi¬ 
cates the typical error on the SINS points. The lensed arcs, as well 
as the WiggleZ and SHiZELS samples, are similar to z ~ 1 main 
sequence star-forming galaxies, though the lensed arcs probe the 
lower-mass end of the relation. The SINS AO sample covers higher 
specific star formation rates, similar to the z ~ 2.5 star-forming 
main sequence. 


position angle and inclination is discussed in Section 13T1) . 
The result is a series of concentric ellipses with the galacto- 
centric radius in pixel i,j given by 


x sin — y cos ^ /rA 
^0 ) ’ U 

where x — \i — i c \ and y — \j—j c \- We then convolve this 
array with the flux-weighted source plane PSF, to account 
for the fact that the preferential direction of magnification 
causes the source plane images to appear smeared in one 
direction. 

Starting from the dynamical centre, we then sum the 
flux within contours of constant r, incrementing r until half 
the total flux is enclosed. This then gives the half-light radius 
7*1/2 given in Table [U The half-light sizes we find are in 
the range ri/ 2 ~ 0.2 — 2.1 kpc. Where we have rest-frame 
optical imaging of a galaxy we carry out the same procedure 
and find that the results are co nsistent with in the quoted 
errors. The half-light radii of the I Jones et al.l (2010) sample 
are r*i/ 2 = 0.5 — 3.5 kpc. As noted in Section [2.41 we find 
galaxies that are systematically smaller than those observed 
in unlensed samples due to the lensing magnification. 


r‘i,j = \ (x cos <f> + y sin 0) + 


3 RESULTS AND ANALYSIS 

3.1 Dynamics 

3.1.1 Dynamical properties of the sample 

From our best estimates of kinematic classification given 
in Table [3] we estimate that 10/17 (59%) of the combined 
lensed sample are rotating discs, 5/17 (29%) are probable 
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Figure 4. Left: Relation between rotational velocity and half-light radius. Right: v ro t/a against half-light radius. The lensed arcs are 
colour-coded by redshift on the same scale in both plots, with filled symbols representing those galaxies identified as dis cs, an d unfilled 
symbols mergers or undetermined syste ms. Comparison samp les are inclu ded from the high -z le nsed sources of lJones et al.l (hold) . and the 
unlensed high-z IFU surveys SHiZELS (ISwinbank et al.ll2Q12h . WiggleZ (IWisnioski et alJl201 lh and the SINS AO sample (Newmar^^ah 
|2(^|). The grey crosses indicate the average error on the SINS points. The dashed line is the fit to the SINS data from lNewman et al. 
( 2013i l. The intrinsically smaller lensed galaxies extend the non-lensed samples to smaller sizes, but the correlation is weaker. 


mergers, and the remaining two are undetermined. The 
merger fraction in our sample is very close to those in un¬ 
lensed surveys, typically 1/3, but we find a higher frac¬ 
tion of rotating s ystems , w hich ty pically make u p 33-44% 
of ot h er samples dForster Schreiber et al.ll2009l : lLaw et all 
120091: IWrieht et alJ 120091 : lEpinat et alJ 1201217 While the 
merger classification adopted is relatively insensitive to spa¬ 
tial resolution, it is likely that many of the galaxies classed 
as ‘dispersion-dominated’ in unlensed studies would present 
with velo city gradients g i ven h igher spatial resolution. For 
example, Ne wman et all (|2013h found that when the same 
34 galaxies observed in seeing-limited conditions were re¬ 
observed with adaptive optics, the fraction classified as 
dispersion-dominated fell from 41% to 6 — 9%. The two 
galaxies we classify as undetermined are A1835-arc7.1 and 
MS 1621-system 1. The former lies on the critical line so it is 
possible that we have not observed the whole velocity gradi¬ 
ent in this source. Alternatively, the turnover in the rotation 
curves of these galaxies might be missed due to the presence 
of a low surface brightness disc. 

The quantity v TO t /cr is a measure of the rotational sup¬ 
port in a system. Galaxies exist on a continuum in v rot /a ; 
they all have some degree of rotation, but those with low 
I’rot /cr values have a larger fraction of their kinematic sup¬ 
port from random motion rather than ordered rotation. Gen¬ 
erally, a simple cut is applied, wher e galaxies with v ro t/cr > 
0.4 are considered rotating systems ([Forster Schreiber et al.l 
2009). Of our combined sample, only A1835-arc7.1 falls be¬ 
neath this cut, with v ro t/cr = 0.2 =L 0.2. However, as dis¬ 
cussed above, we may not have observed the full veloc¬ 
ity gradient in this galaxy. Three of the galaxies that ap¬ 
pear to be interactions or mergers - A1413-arc2.1a, A1413- 
arc2.1b and RXJ1720-arc 1.1+1.2 - all have marginal v ro t/a 
within la of the boundary between rotation-dominated 


and dispersion-dominated systems. The remaining possible 
merger - MS 1621 - has high v ro t/cr = 3 + 1, which could 
be due to a velocity difference between two merging compo¬ 
nents. 

In non-len s ed ga laxies observed with adaptive optics, 
Newman e t al • (2013) found a correlation between v rot /a 
and the half-light radius tt/ 2 . This was suggested to be due 
to a relationship between v ro t and r*i/ 2 , with larger galax¬ 
ies having higher rotation velocities. We add our sample to 
these relationships in Figure 2] As gravitational lensing al¬ 
lows us to observe intrinsically smaller galaxies, we are able 
to extend these relations to smaller sizes, and we find much 
more scatter than in the unlensed systems, with a Spearman 
rank correlation coefficient of p — 0.61 (where p — 0 indi¬ 
cates no correlation, and p — 1 implies perfect correlation). 

Comparing v ro t/a to r*i/ 2 in Figure S] gives similar re¬ 
sults; we find higher scatter at low ?T/ 2 , but there remains a 
weak positive correlation (Spearman rank correlation coef¬ 
ficient p — 0.52) in the combined lensed and unlensed data. 

In both cases, there is no systematic difference in this 
scatter between the rotation-dominated galaxies and those 
identified as non-disc systems. The lack of correlation seen 
in the lensed sample with respect to the unlensed galax¬ 
ies therefore implies a fundamental difference in the kine¬ 
matic properties. As discussed below, this is likely due to the 
systematically smaller sizes and lower masses of the lensed 
galaxies. 

In Figure [5] we plot stellar mass, M*, as a function 
of Vrot/cr, finding that our combined sample continues the 
relation observed in the non-lensed SINS+AO sample to 
lower stellar masses, although with significant scatter. How¬ 
ever, there is little correlation with the unlensed SHiZELS 
and WiggleZ data. This could be because they are observed 
with similar resolution to the SINS data, but have system- 












































































Resolved spectroscopy of gravitationally lensed galaxies 17 




Figure 5. Relation between stellar mass, M*, and v ro t/cr. The 
rotational velocity u ro t is corrected for inclination, and a is the 
luminosity-weighted mean veloci ty dispersion, correcte d for beam 
smearing. The SINS A O sample jNewma n et al.||2013l), SHiZ ELS 
(|Swinbank et al.ll2012l ) and WiggleZ (IWisnioski et alil201lh are 
also shown for comparison with higher stellar mass systems at 
high redshift. The grey cross indicates the typical error on the 
SINS data, and the dashed line is the best fit to the SINS + lensed 
samples. We find that lensed galaxies with high v ro t /& tend to 
have high M*, indicating that galaxies become more dynamically 
‘ settled’ as they build up more mass. 


atically smaller sizes, and thus a larger contribution from 
beam smearing. The best fit to the combined lensed and 
SINS+AO unlensed data is log (M*) = 1.2 log (v TO t/cr) + 9.7. 

If we interpret v ro t/cr as a measure of the ‘order’ in 
a system, where v TO t measures the ordered rotation and a 
constitutes turbulent or disordered motion, this relationship 
suggests that galaxies become more ordered as they build up 
higher stellar masse s, ajorocess described as ‘kinematic set¬ 
tling’ (|Kassin et alJ 12012 ). Therefore, in the lensed sample, 
which probes systematically smaller intrinsic sizes and lower 
masses than the unlensed surveys, we see less contribution 
from ordered rotation as shown in Figure [H and hence do 
not see the correlations exhibited by the larger, more mas¬ 
sive unlensed galaxies. 


3.1.2 The stellar mass Tully-Fisher relation 

The Tully-Fisher relation relates the stellar content of galax¬ 
ies to their rotational velocity (|Tullv &; Fished Il977 ). In a 
model in which gas cools from a dark matter halo into a ro¬ 
tating disc, maintaining the angular momentum of its parent 
halo, the Tully-Fisher relation is interpreted as a relation¬ 
ship between the baryonic content of galaxies and the angu¬ 
lar momentum of their dark matter halos. As such, it is a key 
parameter that models of galaxy evolution must reproduce. 

Attempts to place observational constraints on the red- 
shift evolution of the Tully-Fisher relation have had incon¬ 
clusive results. Observations of the 5-band Tully-Fisher re¬ 
lation at high redshift demonstrate modest evolution, but 
are affected by recent star formation and so evolution of 


Figure 6. The stellar mass Tully-Fisher relation. The stel¬ 
lar mass, M*, is estimated from SED fitting to the available 
HST and Spitzer imaging. The rotation velocity V 2.2 is mea¬ 
sured from the model rotation curves interpolated at a radius 
of r = 2.2 r 1 / 9 , where f jj 2 is the half-light radius. The z = 0 
data are from lPizagno et al.l (|2005l) . and the solid line is the best 
fit to these points. H igh-redshift comparison data are from the 
z ~ 2 lensed arcs oflJones et al.l (120101 ) , the z = 1 lensed arcs 
of ISwinbank et al.l (l2006|), th e DEI MOS z ~ 0.6 and z ~ 1.3 
samples of iMille r e^ak|j^|^ylj2^2 ), the z ~ 2 — 3.5 S INS and 
AMAZE surveys Jcresei et al 120091 ? Gnerucci et al.l201lh a nd the 
SHiZELS 0.84 < z < 2.3 sample of ISwinbank et al.l <20l3 ). The 
dotted line indicates the predicted evolutio n of the Tully-Fisher 
relation at z = 2 from simulations (|Crain et al.ll2009h . We colour- 
code our lensed arcs by redshift, and show rotating disc-like sys¬ 
tems with filled symbols, and mergers and dispersion-dominated 
systems with open symbols. We find that most of our sample is 
largely consistent with the z = 0 relation, but with significant 
scatter discussed further in the text. 


the B-band mass-to-light ratio (|Vogt et al.lll996l . Il997l ). We 
therefore concentrate on the more physical stellar mass 
Tully-Fisher relation, but this too has produced mixed re- 


dConselice et al. 2005; [Flores et al. 2006; Cresci et al. 

20061; 

Kassin et al.l 2007 

; Puech et al. 20081: IVergani et al. 

2012; 

Miller et al. 201 ll. 

20121; Swinbank et al. 20121). Miller et al. 


that show the greatest offset from the local relation, which 
they interpret as evidence that these galaxes have yet to 
‘mature’ onto the Tully-Fisher relation, per haps related to 
highe r gas fractions in their discs. Similarly, I Gnerucci et al.l 
ll201lll investigated the stellar mass Tully-Fisher relation at 
z ~ 3 and found marginal evolution but with a large amount 
of scatter, concluding that the relation may not be in place 
at this epoch. 

To compare different galaxies, it is important to mea¬ 
sure r’rot in a consistent manner. Measurements of the rota¬ 
tion curve are naturally sensitive to the surface bright ness 
of the galaxy and the depth of the observations. ICourteaul 
( 1997 ) compared various methods of measuring v ro t and 
found that V 2.2 (the velocity interpolated from the model 
rotation curve at a radius r m 2.2 h, where h is the disc scale 
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Figure 7. The ratio between stellar mass and galaxy radius as a 
function of rotational velocity, where radius is defined as 2.2 disc 
scale lengths and the velocity is corrected for inclination. Compar¬ 
ison data is from the same sources as in Figure [6] The dashed line 
indicates v 2 = GM*/r. The lensed galaxies and some data from 
the other samples, in particular from SHiZELS, lie along this line, 
indicating that their dynamics are dominated by baryons in the 
discs. There is no apparent systematic difference between disc¬ 
like (filled symbols) and mergers/ dispersion-dominated systems 
(open symbols) in the lensed samples. Offsets from this line are 
likely related to the gas fractions of the galaxies. 

length) performed best in terms of minimal internal scatter 
and residuals from the Tully-Fisher relation and provided 
the best match to radio (21cm) results. In high-redshift 
galaxies with clumpy and disturbed morphologies, the light 
does not follow a pure exponential disc profile, so we can¬ 
not strictly use this definition with our data. We therefore 
adopt the half-light radius, 7T/ 2 , as a proxy for the disc scale 
length, and measure V2.2 from the model rotation curve fits, 
corrected for inclination, at r — 2.2 r 1 / 2 . For consistency 
across the sample, we re-analyse the velocity fields of the 
I Jones et al.l |2010 ) sample in order to extract V2.2 from these 
galaxies. Our derived V2.2 values are slightly lower than the 
inclination-corrected rwx that they derive by an average of 
20 %. 

We plot the stellar mass Tully-Fisher relation in Fig¬ 
ure [6] We find that most of our lensed arcs are consistent 
with the Tully-Fisher relation observed in other samples, 
but with considerable scatter. Three galaxies he well above 
the local relation, with low V2.2 for their stellar masses. Of 
these, two are from the I Jones et ID 1 I 2010 I ) sample; these are 
C10949-arcl and the Cosmic Ey e. The former was identi¬ 
fied as a merger by Ijones et~al ] fcoinh . though th e latter 
appe ars to be a rotationally-supported disc (see lStark et al.l 
120088 ). The remaining galaxy lying above the local relation 
A1413-arc2.1a, which appears to be interacting with A1413- 
arc2.1b. 

Simulations predict that the zero-point of the Tully- 
Fisher relation should decrease at high redshift (e.g. 
iMcCarthv et alJl2012l : ICrain et al.ll200d indicated by a dot¬ 
ted line in Figure[6]). Semi-analytic models in which the disc 
dynamics are driven by the dark matter halo make a simi¬ 


Figure 8. Relation between Star formation rate (SFR) in 
clumps and their sizes. The z = 0 sample from the SINGS 
survey (|Kemncutt__eL_al] [20 03|) a nd the z ~ 1 — 1.5 lensed 
arcs from iLivermor^e^il .1 J201 2 ah are extracted from Ha nar¬ 
rowband imaging. W e als o show the z y 5 le nsed arcs of 
ISwinbank et al.l (|2007l ) and ISwinban k e^ak|_j[2009|)_ and th e un- 
lensed hig h -z samples from SHiZE LS dSwinbank et al.ll2012h and 
WiggleZ (IWisnioski et al.l l2012al j for comparison. The clumps 
from our sample of lensed arcs are colour-coded by redshift. The 
dashed lines are contours of constant surface brightness with the 
zero-points fit in four redshift bins. We find evolution in this zero- 
point such that high-redshift galaxies host clumps of higher sur¬ 
face brightnesses. 

lar prediction ([Dutton et al.ll201lh . However. iBenso D (l2012h 
predicts a mo dest positi v e incr ease in the zero-point at high 
redshift. As lMiller et al.l ( 20121 ) point out, the primary differ¬ 
ence is that the lBensonl model has a larger contribution from 
baryons. High-redshift galaxies are intrinsically smaller than 
their present-day counterparts - especially in the lensed sam¬ 
ple that probes the smaller, fainter population - and so we 
might expect the d ynamics with in 2.2r-| / 2 to be dominated 
by baryons (see also lMiller et ahlfeoilh . In the simplest case 
of orbital circular motion, we should see v 2 = GM/r for 
total enclosed mass M. To test whether this is a better de¬ 
scription of our data, we divide the y-axis of Figure [6] by 
2.2 t*i/ 2 and show the result in Figure 0 The scatter is sig¬ 
nificantly reduced in comparison to Figure [6j and the line 
v 2 = GM*/r provides a good fit to the lensed and SHiZELS 
data, although the same three lensed galaxies remain far 
outside this relation. The varying gas fractions in the galax¬ 
ies will contribute significantly to the scatter in this relation. 
The other comparison samples have a systematically higher 
velocities at a given M*/r, perhaps suggesting that their 
dynamics do have contributions from the dark matter halo, 
whereas the dynamics of the lensed galaxies within the cen¬ 
tral few kpc are dominated by baryons. 

3.2 Star-forming clumps 

With the high spatial resolution afforded by gravita¬ 
tional lensing, it is possible to extract star-forming 
clumps from lensed galaxies on scales comparable to 
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Table 5. Properties of the star-forming clumps 


Name 

Notes: 

radius 

pc 

(a) 

SFR 
M@yr —1 

(b) 

a 

km s _1 
( c ) 

M ACS0744-system3-1 

300 ± 100 

0.26 ± 0.02 

70 ± 7 

MACS0744-system3-2 

300 ± 100 

0.18 d= 0.01 

84 d= 8 

MACS0744-system3-3 

160 d= 30 

0.04 d= 0.01 

51 d= 8 

MACS 1149-arc A 1.1-1 

450 ± 70 

0.51 d= 0.05 

74 d= 10 

MACS1149-arcAl.l-2 

150 d= 50 

0.15 d= 0.02 

120 d= 19 

MACS1149-arcAl.l-3 

160 d= 30 

0.06 =b 0.01 

134 d= 31 

MACS1149-arcAl.l-4 

210 d= 50 

0.12 ± 0.02 

80 ± 17 

MACS1149-arcAl.l-5 

180 ± 70 

0.07 d= 0.01 

95 d= 39 

MACS1149-arcAl.l-6 

150 d= 20 

0.05 d= 0.01 

71 d= 17 

MACS 1149-arcAl. 1-7 

160 d= 40 

0.04 ± 0.01 

123 d= 33 

MACS1149-arcAl.l-8 

140 d= 50 

0.03 d= 0.01 

115 ± 30 

MACS1149-arcAl.l-9 

140 ± 20 

0.02 d= 0.01 

122 d= 61 

MACS0451-system7-l 

500 ± 200 

1.2 d= 0.4 

87 d= 6 

MACS0451-system7-2 

300 ± 200 

0.5 ± 0.2 

84 d= 6 

A1413-arc2.1a-l 

700 ± 300 

2.3 d= 0.2 

53 d= 8 

A1413-arc2.1b-l 

200 ± 100 

0.18 d= 0.02 

57 d= 9 

A1413-arc2.1b-2 

260 d= 60 

0.15 d= 0.02 

58 d= 10 

A1413-arc2.1b-3 

300 ± 100 

0.26 d= 0.03 

54 d= 9 

A1413-arc2. lb-4 

140 ± 40 

0.02 d= 0.01 

84 ± 32 

A1835-arc7.1-l 

90 d= 60 

0.10 d= 0.01 

77 d= 14 

A1835-arc7.1-2 

200 ± 100 

0.12 d= 0.02 

67 d= 13 

MS1621-systeml-l 

900 ± 400 

3.1 ± 0.3 

63 d= 9 

MS1621-systeml-2 

600 d= 200 

0.8 d= 0.1 

66 d= 11 

MS1621-systeml-3 

340 d= 70 

0.18 d= 0.03 

55 d= 12 

MS1621-systeml-4 

200 d= 100 

0.08 d= 0.01 

58 d= 12 

RXJ1720-arcl.l+1.2-l 

700 ± 200 

0.77 d= 0.08 

61 d= 8 

RXJ1720-arcl.l+1.2-2 

600 ± 100 

0.59 d= 0.06 

67 d= 10 

RXJ1720-arcl.l+1.2-3 

400 ± 100 

0.28 ± 0.03 

76 ± 11 

RXJ1720-arcl.l+1.2-4 

600 ± 100 

0.65 d= 0.06 

65 d= 9 

RXJ1720-arcl.l+1.2-5 

300 ± 100 

0.23 d= 0.03 

67 d= 12 

A1689-arc2.1-l 

180 d= 60 

0.2 d= 0.1 

116 d= 45 

A1689-arcl.2-l 

150 d= 80 

0.13 d= 0.04 

103 d= 16 

A2895-arcl.l+1.2-l 

400 ± 100 

1.8 ± 0.7 

63 d= 4 

A2895-arcl.l+1.2-2 

500 ± 300 

2.3 ± 0.8 

61 d= 4 

A2895-arcl.l+1.2-3 

60 d= 10 

0.03 d= 0.01 

57 d= 5 

A2895-arc2.2-l 

290 d= 90 

5. d= 2. 

51 d= 3 

C10024-arcl.l-l 

420 ± 70 

5. d= 1. 

65 ± 11 

C10024-arcl.l-2 

600 ± 200 

9. d= 2. 

60 d= 9 

C10024-arcl.l-3 

370 d= 50 

4. d= 1. 

63 d= 14 

C10024-arcl.l-4 

400 d= 100 

3.4 d= 0.8 

62 d= 13 

C10949-arcl-l 

1000 ± 100 

16. ± 4. 

72 d= 13 

C10949-arcl-2 

800 ± 300 

10. d= 2. 

63 d= 8 

C10949-arcl-3 

800 ± 300 

11. d= 3. 

56 d= 6 

C10949-arcl-4 

700 ± 100 

8. ± 2. 

76 d= 14 

M ACS0712-systeml-1 

220 ± 70 

3.3 ± 0.8 

74 d= 11 

MACSJ0744-arcl-l 

300 ± 100 

0.7 d= 0.2 

93 d= 26 

MACSJ0744-arcl-2 

300 ± 100 

0.4 ± 0.1 

136 d= 45 

MACSJ0744-arcl-3 

300 d= 100 

0.5 d= 0.1 

118 d= 14 

CosmicEye-1 

600 ± 100 

25. ± 6. 

41 d= 15 

CosmicEye-2 

500 d= 100 

12. d= 3. 

40 ± 14 


Notes: All quantities are corrected for lensing where appropriate. 
(a): The clump radius is calculated as described in the text, using 
the area enclosed by an isophote and assuming circular symmetry. 
The error bar encompasses an alternative definition using the 
FWHM of the clump emission line profile. (6): SFR is calculated 
from the total Ha or H/3 flux enclosed in the isophote and is 
corrected for extinction, (c): The velocity dispersion, cr, is the 
mean local value within the clump and is deconvolved for the 
local velocity gradient. 


giant Hu regions in local galaxies. It has long been 
noted that galaxies in the high-redshift Universe tend 
to be ‘clumpier’ then their local counterparts, with 
clumps often dominating the host galaxies’ morpholo¬ 
gies and givi ng rise to ‘clump cluster’ and ‘ chain ’ 
galaxies (e.g. Cowie et al. Il995l : lElmegreen et al.l [2004; 
lElmegreen fe Elmegreenl l2005h . Numerous high-redshift 
galaxy surveys have isolated clumps jmd analysed their 
properties on k p c scales dForster Schreiber et al.l l201ll 


Genz el et al.l l201ll : IWisnioski et al.l 2012al : ISwinbank et al l 


2012), but to compare them directly to star formation in 


the local Universe, which takes place in Hu regions on scales 
of ~ 100 pc, the sp atial resolution afforded by gravitational 
lensin g is required (ISwinbank et al.ll201 ll . 120091 : 1 Jones et al.l 


lensi ng is required (bwmbank < 
120101 : [Livermore et al.l 2012ah . 


Different studies of star-forming clumps have used 
a variety of methods to detect and isolate them from 
their host galaxies. Discussion s of the meri t s of va rious 
techniques can b e foun d in IWisnioski et al.l (|2012al i and 
ILivermore et al. (l2012alh W e adopt the same method as 


Livermore et al.f 


■ , iip— ii .uiTMiii 1 iim ' an d use the 2D version of 

CLUMPFIND ([Williams et alJ Il994h . This routine uses mul¬ 
tiple isophotes, starting by defining clumps in the brightest 
regions and then moving down through the isophote levels. 
Any isolated contours are defined as new clumps, and any 
which enclose an existing peak are allocated to that clump. 
A contour which encloses two or more existing peaks has 
its pixels divided between them using a ‘friends-of-friends’ 
algorithm. We define the lowest contour at 3 < 7 , where a is 
the standard deviation of pixels in the underlying galaxy. 
We then add additional contours in la intervals up to the 
peak in the image. _ _ 

For IFU data, iGenzel et al.l (2011) advocate defining 
clumps in slices in velocity. We test this method with our 
data and find that it produces similar results to using 
the integrated intensity maps, and where it does produce 
clumps that are not evident in the integrated maps, they 
are marginal detections found only in the lowest (3 — 4a) 
contour levels. We therefore use the integrated Ha (or H/3) 
maps to define clumps. This has the added advantage of be¬ 
ing the method most easily com parable to the narrowband 
imaging in ILivermore et al l ()2012al i, where we used maps of 
Ha integrated over a wide velocity range. 

To ensure self-consistency across the entire sam¬ 
ple of lensed a rc s, we reanalyse the intensity maps of 
the I Jones et ah (2010) samp le to extract clumps with 
CLUMPFIND. Ijones et al.l (|2010h used a single isophote per 
galaxy to extract t heir clump properties, wh ereas the algo¬ 
rithm developed bv ILivermore et ahl (|2012aT ) used the noise 
properties of the image to robustly define a series of isophote 
levels in a way that does not require any fine tuning. The 
variation in methods means that the clumps extracted are 
not identical; a lower isophote will increase both the size 
and luminosity of the extracted clump, and vice versa. Gen - 
erally, we use lower isophote levels than I Jones et al ] LoioL 
leading to larger sizes (but within ~ la of the published 
values) with correspondingly higher luminosities. The pri¬ 
mary difference is that we also extract some smaller, lower- 
luminosity clumps that are missed when adopting a single 
isophote. 

Once the clumps are defined, we sum the spectra in their 
component pixels and fit a Gaussian emission line (a single 




































































20 R. C. Livermore et al. 


line for H a or a triplet for H/3 and [Om]) to obtain their 
intensity and velocity dispersion. For those galaxies observed 
in H/3, we estimate the equi valent Ho intensity b y assuming 
case B recombination and a lCalzetti et al.l ( 2000 ) extinction 
curve based on the attenuation Ay estimated from the SED 
fitting. No additional extinction is assumed in the nebular 
lines. Thus, the estimated Ho luminosity, L Ha is given by 

•f Ha = 2.86 X 10 7 - 96 Lh/3* (3) 

From Lh oi , we ca lculate the SF R using the presc ription 
of iKennicuttl (|l998al ) adapted to a IChabriei] (2003) initial 
mass function (IMF) and extinction Aua = 0.82 Av- 

The effective radius, r, of the clump is estimated by 
summing the area of all pixels contained in the clump, A. 
We then subtract the area of the effective source plane PSF 
in quadrature and define an effective radius 


(4) 


We identify 50 clumps in our sample, with between 1 
and 9 clumps per galaxy, and list the derived properties in 
Table [5] 

I Jones et al ] fcoich found that the clumps in high- 
redshift lensed galaxies had surface brightnesses that are 
higher than those of local Hu regions by up to 100 x. Such 
intense star-forming regions a re found locally in merging 
systems such as the Antennae ([Bastian et al.ll2006h , but ap¬ 
peared to be ubiquitious in isol ated discs at high-z. Simi- 
lar results were found at z ~ 5 (jSwinbank et al. 20091) and 
in kpc - scale clumps from unle nsed sources ( Wisnioski et al.l 
l2012bl : ISwinbank et al.l l2012h . To invest igate the o rigin 
of thes e high-surface brightness clumps, iLivermore et al.l 
(|2012al ) observed 57 clumps in 8 galaxies with H a narrow- 
band imaging, and demonstrated that their surface bright¬ 
nesses evolve with redshift. In Figure [8j we expand on 
this relation adding our new clumps. We note that one 
galaxy - MACS 1149-arcA 1.1 - is included in both the current 
IFU sample and the narrowband imaging sample. Only one 
clump (the galaxy bulge) was observed in the narrowband 
sample, and we confirm that its luminosity is consistent with 
the previous observation within the measurement errors. 

With the new data, the evolution in clump surface 
brightnesses is confirmed and extended out to z > 3. We 
show the mean surface brightness, U c i ump , in four redshift 
bins in Figured] The best fit to the redshift evolution of the 
clump surface brightness is 


log ( - Sd "T _ 2 ) = (3.5 ± 0.5) log (1 + *)-( 1.7 ± 0.2). 

\M© yr -1 kpc V 

(5) 

Clearly there is a selection effect involved, as we would 
not observe low-surface-brightness clumps in the high -z 
galaxies. However, there appears to be evolution in the prop¬ 
erties of the brightest clumps, as these are not seen in iso¬ 
lated local galaxies but seem to be ubiquitous at high red¬ 
shift. We will discuss this further in the context of the clump 
luminosity function along with implications of this evolution 
in Section PI 

The IFU data of our current sample allows us to add the 


extra dimension of velocity dispersion, cr, to this analysis. If 
the clumps are gravitationally bound , they sh ould f ollow a 
relation of the form L oc cr 4 ([Terlevich fc Melnickl ll98~ l). 
However, not all studi es of local star-fo rming regions have 
found this relation; lArsenault et al.l ([19901 ) found that it ap¬ 
plied only to the brightest, ‘first-ranked’ regions in galaxies, 
with no L — a or r — a relation applying to their entire sam¬ 
ple. Others have found a shallower L — a relation indicative 
of density-bounded regions, where only a fraction of the ion - 
ising photons are able to escape ([Rozas et al.l l200ll . 120061 ) . 
It is not clear that clumps should be expected to be viri- 
alised; indeed, Giant Molecular Clouds (GMCs) lo cally are 
not necessarily gravitationally bound (iDobbs et al.ll201 lh . 

In Figure [9] we plot a as a function of both the clump 
size and lumi nosity. As a guid e, we overlay the relations 
found bv ISwinbank et all (2012 ): 

108 ( kjs) = < 101 ± »• “>« (i^r) + <°' 8 ± °» fe > 
los (s£r) = < 8 ' 81 ± (l^r) + < 34 ' 7 ± °' 4 > • 

(7) 

Although our data scatter around these relations, we do 
not find any correlation between a in our clumps and either 
r or L. 

We first consider whether this is an effect of spectral 
resolution. We deconvolve a for the instrumental resolution 
as measured from sky emission lines, but there is some lower 
limit at which we will not be able to measure the broadening 
of the line. To test where this limit is, we construct a set 
of 1000 Gaussian emission lines with varying widths of < 
100 km s -1 and add noise so that the final signal-to-noise, 
S/N= 5 (this gives us a conservative estimate, as S/N> 
5 is the constraint we set to detect an emission line). We 
then convolve the resulting spectrum with the instrumental 
resolution, and apply the same emission line fitting routine 
used on the data to recover the line properties. We find that 
we consistently recover the input line width to within 20% 
at cr > 53kms _1 . We show this line on Figure [9] 

The clumps that lie close to this line may be affected 
by the resolution limit; however, there remains an excess of 
clumps at high-cr that are well-resolved and lie above the ex¬ 
pected relations. It therefore appears that the clumps are not 
necessarily virialised, and the observed a may contain con¬ 
tributions from a range of processes including gravitational 
instabilities and star formation feedback. A further possible 
contribution to the high cr values could come from the super¬ 
position of multiple clumps along the line of sight. This is an 
additional source of uncertainty in all values derived from 
clumps, and motivates studies with higher spectral resolu¬ 
tion so that any such multiple clumps can be kinematically 
separated. 


4 DISCUSSION 
4.1 Clump formation 

I Livermore et al.l (|2012al ) demonstrated that an effective 
means of quantifying the ‘clumpiness’ of a galaxy is the 
clump luminosity function. The luminosity function can 
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in the text). The addition of lensed galaxies extends the high -z samples to smaller sizes and lower luminosities while cr is similar to the 
unlensed samples. Although partially an effect of the spectral resolution limit, we find an excess of clumps with high resolved cr for their 
sizes and luminosities. 


be described by a Schechter function with a cut-off at 
some maximum luminosity Lo, which describes the bright¬ 
est clumps. The presence of more bright clumps in a galaxy 
causes it to appear as ‘clumpy.’ I Livermore et al.1 ( 2012ah 
showed that Lo increases with redshift, giving rise to the 
appearance of more clumpy galaxies at high-redshift. 

We show the luminosity functions in Figure [10] with the 
new data added. We account for the varying surface bright¬ 
ness limits of the observations by normalising each bin by 
the number of galaxies in which we should be able to identify 
clumps of that luminosity, and the error bars represent the 
Poisson error from counting clumps. The z < 1.5 bin com¬ 
bines the two z < 1.5 galaxi es from the current sam ple with 
the seven galaxies from the lLivermore et al.l (|2012al ) narrow- 
band imaging sample (MACS 1149-arcA 1.1 is removed from 
the latter as it is duplicated in the current sample), and 
the z > 3 bin i nclud es th e two lensed z y 5 g alaxies of 
ISwinbank et al.l (|2007h and lSwinbank et al.l (|2009T ). 

To each redshift bin, we fit a Schechter function of the 
form 


w(>i >= Ar °(0’ < “ p (iF (8) 

w here the powe r-law slope is fixed to a = —0.75 from 
iHookins et al.l ( 2012 ). We do not attempt to fit the faint- 
end slope as the sharp turnovers evident in Figure [TO] are 
likely due to incompleteness and small number statistics in 
the faintest clumps. The normalisation No is arbitrary, so 
we fit No to the z = 0 data. We then keep No and a fixed 


while allowing the cut-off Lo to vary, and we find the best-fit 
Lo for each bin with a % 2 minimisation procedure. The best 
fits are 


i 2 *5±l;l x 10 41 , z = 0 
4.0i°;^xl0 41 , 1 < z < 1.5 

4+% x 10 42 , 1.5 < z < 3 

1.3tg;jj x 10 43 , z > 3 


where the errors are estimated using a bootstrap 
method. 

We therefore find that the cut-off of the clump lumi¬ 
nosity function evolves with redshift, such that high-redshift 
galaxies have more high-luminosity clumps. The best fit to 
the redshift evolution of L 0 is 


log -~ t = (2.0 =b 0.7) log (1 + z) + (41.0 =b 0.2) . (9) 

\ergs- 1 / 

In combination with the evolution in clump surface 
brightness discussed in the previous section, it follows that 
high-redshift galaxies tend to have a higher number of 
bright, high surface-brightness star-forming clumps, which 
come to dominate the galaxy morphology and therefore give 
rise to the population of ‘clumpy’ galaxies seen at high red- 
shifty _ 

I Livermore et al.l (|2012aT ) suggest that the cause of these 
bright clumps is a combination of high gas fractions at high 
redshift and the evolving galaxy dynamics. Briefly, the ar¬ 
gument is based on the assumption that clumps form in 
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Figure 10. Cumulative clump lu minosity functions. The z = 0 
data are from the SING S survey (| Kennic utt et al.ll2003l) , with 
clumps extracted by ILivermore et al. I (l2012al) . The z < 1.5 
bin inc ludes the narrowband imaging sample of ILivermore et al.1 
(2012a) combined with the two z < 1.5 galaxies from this work. 
Each luminosity bin is normalised by the number of galaxies con¬ 
tributing to that bin, accounting for variations in the depth of ob¬ 
servations as described in the text. The dashed lines are Schechter 
f unction fits, where the z = 0 fit is from the simulations of 
[Hopkins et al ] (2012), and we keep the same normalisation but 
allow the cutoff to vary to fit the high-redshift data. We find that 
this cut-off evolves with redshift. 


marginally stable discs described by a iToomrel (1 19641 ) pa¬ 
rameter Q where 


<io) 

for epicyclic freque ncy n » a/2 v/ r and velocity v at 
galactocentric radius r. iHopkins et al.1 ( 20121 ) argued that 
star-forming galaxies tend to self-regulate to maintain Q ~ 
1; instabilities cause the gas to collapse to form stars, while 
feedback from the star formation stabilises the disc. 

The mass required for collapse on a scale R is the Jeans’ 
mass, Mj, given by 


Mj = (11) 

If Mj is related to the maximum mass of clumps in 
a galaxy, Mp (i . e. the cut-off of the mass function), then 
ILivermore et al.l (|2012al ) showed that combining Equations 
cm and m with an assumed turbulent power spectrum 
gives 


Mo - 


37r 3 G 2 Eg 

2 K 4 ' 


( 12 ) 


Thus, Mo depends on the disc surface density and 
epicyclic frequency. To relate this predicted mass to the 
observable Lh«, we use an empirical relation between 


the Ha-derived SFR and mass of local molecular clouds, 
SFR(M 0 yr“ 1 ) = 4.6 ± 2.6 x 1O“ S M 0 llLada et, al.l 1 201 (I) . 
This relation is valid for high-density gas, and thus ap¬ 
propriate for star-forming clumps, and is consistent with 
the independently derived SFRs and gas masses of star¬ 
forming clumps reported in a lensed z = 2.3 galaxy by 
ISwinbank et al.l (|201ll ). Adopting this conversion results in 
cut-off masses, Mo, from our clump luminosity functions of 


Mp 

M 0 


'4±* x 10 7 , z = 0 

l.OtJ.a x 10 8 , 1 < z < 1.5 

9t£xl0 8 , 1.5 < z < 3 

^xlO 9 , z > 3. 


We can also relate the SFR surface density, T,sfr , of 
clumps to the disc surface density, Eo, using the Kennicutt- 
Schmidt law. The relationship be tween the surface density 
of gas, E gas , and Esfr found bv iKennicu tt] (Il998bl) is 


Esfr = A ( — V 

Mq yr -1 kpc -2 \M & pc~ 2 J 


(13) 


where A = (2.5 ± 0.7) x 10“ 4 and n = 1.4 ± 0.15. The 
disc surface density Eo is a combination of gas and stars (E g 
and E* respectively^ which contribute differently to the disc 
stability. Following iRafikovI (2001), we use 


£»-£« + (tttj) e - (») 

where f a = cr*/cr g « 2 is the rat io of t he velocity disper- 
son of the stars to that of the gas (jKorchagin et al.l 120031 ). 
We define the gas fraction, / gas , such that 


/gas — 


Mg; 


Eg + E* 


(15) 


if we measure the gas and stars over the same area. 
Combining Equations HU HU and H21 we have 


Esfr = A 



/gas (1 + fcr) 


/gas (1 + fa) 2(1 /gas) 


n 

■ (16) 


Thus, from Equations \T2\ and HU we find that the cut-off 
mass, Mo, and the clump SFR surface density, Esfr, depend 
on Eo, /t and / gas . In order to understand how Mo and Esfr 
evolve with redshift, we therefore need to understand the 
evolution of Eo, k and / gas - 

We can estimate Eo and / gas from our data using the 
Esfr derived from Ha or H/3 and applying the Kennicutt- 
Schmidt law from Equation HU to derive E gas . To estimate 
E*, we use the total M* from the SED fitting and apportion 
the stellar mass according to the fractional flux in each pixel 
of the reddest available HST image. We then combine the 
E gas and E* maps to create a map of Eo. There are many 
uncertainties in estimating of / gas from Ha, specifically, the 
assumption of a constant Kennicutt-Schmidt law and in the 
reliability of the SED-derived stellar masses; nonetheless, we 
note that our estimate for the Cosmic Eye of M gas ~ 10 9 M© 
is very c lose to t he CO - derive d value of M gas = (9.0 ± 1.6) x 
10 8 M© (|Riechers et al.ll2010l ). In our sample, we find a mean 

































































Resolved spectroscopy of gravitationally lensed galaxies 23 


10.00 


o 

1.00 


I 

W° 

0.10 


0.01 

10 7 10 8 10 9 10 10 
M 0 (M 0 ) 

Figure 11. Evolution of clump surface brightness Ssfr and the 
cut-off mass Mo derived from the clumps luminosity functions. 
The dashed line shows a predicted track based on Equations fl2l 
and 1161 assuming evolution in m ass a nd epicyclic frequency from 
the simulations of[Dutton^^ah| |2^n|) and gas fraction evolution 
of /gas oc (1 + z) 2 (|Geach et al.ll201ll ). with f a = 2 . The shaded 
regions indicate the predicted locations of clumps in each redshift 
bin, for a range of gas fractions / gas oc (1 + 2; )( 2 -0±0.5) rj^ 
are consistent with the model, demonstrating that the ‘clumpy’ 
morphologies of high- 2 : galaxies are driven by evolution in their 
gas fractions and dynamics. 
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< /gas >= 0.4 d= 0.1, consistent with the evolution obs erved 
in CO (e.g. iGeach et al.ll201ll : [Livermore et al.ll2012bh . 

We estimate ft by taking the velocity at the half-light 
radius 7 * 1/2 along the dynamical axis from the disc fitting, 
where there is no inclination correction for the radius. This 
measurement is highly uncertain as the measurement of 
rq /2 is affected by the clumpy morphologies and the irreg¬ 
ular source plane PSF, and the measurement of the ve¬ 
locity includes uncertainty in the inclination. We find a 
median ft ~ 90kms _1 , consistent with the predictions of 
iDutton et al. 113) • 

Using these models for the evolution of / gas , Uo and ft, 
we should be able to predict how the star formation surface 
density of clumps, E c i ump and the cut-off mass, Mo, evolve. 
In Figure IIII we present an updated version of the figure in 
I Livermore et al.l (|2012al . Figure 9) including the IFU data. 
We find that the evolution in the cut-off of the clump mass 
function and in the surface brightness of clumps is consis¬ 
tent with the model predictions given by Equations 1121 and 
m given evolution in the galaxy mass and dynamics from 
theoretical predictions, and empirical gas fraction evolution. 
This indicates that galaxies at high redshift appear clumpy 
because their gas-rich discs fragment on larger scales, lead¬ 
ing to star-forming clumps that dominate the morphology of 
the galaxy. Thus, despite the different appearance of high- 
redshift clumpy galaxies to local spirals, we find that they do 
not require a different ‘mode’ of star formation; the clumps 
can arise naturally as a consequence of the larger gas reser¬ 
voirs. 


5 CONCLUSIONS 

We have presented integral held spectroscopy around the Ha 
or H /3 emission lines of 12 gravitationally lensed galaxies at 
1 < * < 4, obtained with VLT/SINFONI, Keck/OSIRIS 
and Gemini/NIFS. We com bine these data with 5 galax¬ 
ies from I Jones et al.l ( 2010 ) and investigate the dynamics 
and star formation properties of 17 high-redshift galaxies. 
The galaxies ah benefit from magnification due to gravita¬ 
tional lensing, increasing the flux by factors of 1.4 — 90 x 
and providing sub-kiloparsec spatial resolution, though we 
stress that uncertainties are introduced by the anisotropic 
magnification. 

Our combined sample has stellar masses M* ~ 6 x 10 8 — 

6 x 10 lo M© and dust extinctions of Ay ~ 0.4 — 1.1. The 
intrinsic star formation rates derived from the Ha (or H /3) 
emission, after correcting for lensing and dust extinction, 
are SFR^ 0.8 — 40M©yr _1 . The use of gravitational lensing 
therefore allows us to probe the representative star-forming 
population at this epoch, with lower M* and SFRs than 
other high- 2 ; surveys. 

We fit rotating disc models to the sample to derive in¬ 
clinations and dynamical axes, from which we extract rota¬ 
tion curves and velocity dispersion profiles. Ah of the galax¬ 
ies in the combined sample have velocity gradients, and 14 
out of 17 are well-fit by an arctan function indicative of 
rotation. Of the sample, 59% have velocity profiles indica¬ 
tive of rotation, centrally-peaked velocity dispersion profiles 
and sufficiently high v/a to be classed as possible rotating 
discs. Of the remainder, 29% have disturbed profiles that 
may be indicative of merging or interacting systems. The 
remaining two galaxies are undetermined as we do not see 
a turnover in their rotation curves. The merger fraction is 
in good agreement with other high- 2 : kinematic surveys (e.g. 


Forster Schreiber et al. 

2009; Law et al. 20091: Wrig 

it et al. 

2009: lEoinat et al. 2012 

L while, like Newman et al. 

(l2013h. 


we find that higher spatial resolution causes a higher frac¬ 
tion of the remainder to be classed as rotating discs. 


As gravitational lensing allows us to probe intrinsically 
smaller galaxies, we extend relationships observed between 
galaxy size, velocity and v/a in unlensed samples to smaller 
sizes. Our results support previous work finding that larger 
galaxies generally have a higher contribution to their kine¬ 
matic support from ordered rotation compared to random 
motions, but we find a large degree of scatter in these smaller 
galaxies. 

We find that the sample is consistent with the local stel¬ 
lar mass Tully-Fisher relation with no coherent evidence fo r 
redshi ft evolution, in common with the work of lMiller et al .1 
ll2012ly We further demonstrate that the rotation in our 
sample, which comprises systematically smaller sizes than 
un lensed studies, could be dominated by baryons (see also 
iMiller et al.ll201lh . 

We extend the w ork of I Jones et 3D (I 2010 I') and 
I Livermore et al.l (|2012al ) by detecting 50 star-forming 
clumps in our sample and study their luminosities, sizes 
and velocity dispersions. In common with previous work, we 
find that the surface brightness evolves with redshift, but 
we extend this evolution to z > 3. The average star forma¬ 
tion density in the brightest clumps evolves with redshift 
as log (E c i ump /M© yr _1 kpc -2 ) = (3.5 =L 0.5) log (1 + z) - 
(1.7 ± 0.2). However, we find that the clumps have similar 
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velocity dispersions to unlensed high -z samples while be¬ 
ing smaller and less luminous; thus, they introduce a much 
larger degree of scatter into the L — a and a — r relations 
observed in other studies. This could be an indication that 
these clumps are not virialised, and that their velocity dis¬ 
persions may have additional contributions from star forma¬ 
tion feedback or gravitational instability. 

We construct luminosity functions of the clumps, and 
find that they can be fit by a Schechter function in which 
the cut-off evolves to higher luminosities at higher red- 
shifts. We find that the cut-off evolves as log (Lo/ergs _1 ) = 
(2.0 db 0.7) log (1 + z) + (41.0 =b 0.2). This supports the pic¬ 
ture in which ‘clumpy’ galaxies arise at high-z because gas- 
rich, turbulent discs fragment on larger scales, resulting in 
star-forming regions large enough to dominate the morphol¬ 
ogy of the galaxy. The gas content of high-redshift galaxies 
remains relatively unexplored due to the small number of 
observations of molecular gas in normal high-redshift star¬ 
forming galaxies. It is to be hoped that the redshift evolu¬ 
tion of gas fractions and the universality of the Kennicutt- 
Schmidt law can be tested as ALMA in full science opera¬ 
tions opens up the more normal star-forming population. 
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Table Al. HST optical photometry used for SED fitting. All fluxes are given in njy and are corrected for lensing. Upper limits are 
given at the 3 cr level. 


Name 

F225VE 

F275VE 

F336VE 

F390W 

F435VE 

F475VE 

F555VE 

F606W 

F625VE 

F775W 

MACS0744-system3 

50 ±10 

80 ±20 

250 ± 50 

160 ± 30 


170 ± 30 

200 ± 40 

230 ± 40 


410 ± 80 

MACS1149-arcAl.l 

180 ± 30 

240 ± 50 

380 ± 70 

370 ± 70 

450 ± 90 

430 ± 80 

390 ± 70 

420 ± 80 

450 ± 80 

440 ± 80 

MACS0451-system7 

A1413-arc2.1a 








150 ± 30 


300 ± 60 

A1689-arc2.1 


0.7 ±0.1 

2.7 ±0.5 



14 ±3 



22 ±4 

29 ±6 

A1689-arcl.2 


0.39 ± 0.07 




5 ± 1 



11 ±2 

11 ±2 


Table A2. HST near-infrared photometry used for SED fitting. All fluxes are given in njy and are corrected for lensing. Upper limits 
are given at the 3 a level. 


Name 

F814VE 

F850LP 

F105VE 

F110W 

F125VE 

FU0W 

F160W 

MACS0744-system3 

450 ± 80 


700 ± 100 

700 ± 100 

480 ± 90 



MACS 1149-arcAl. 1 

500 ± 100 

700 ± 100 

810 ± 150 

880 ± 160 

870 ± 160 

920 ± 170 

1000 ± 200 

MACS0451-system7 

A1413-arc2.1a 

160 ± 30 

360 ± 70 


280 ± 50 


450 ± 80 


A1689-arc2.1 


40 ±8 



300 ± 60 


420 ± 80 

A1689-arcl.2 


13 ±2 



100 ± 20 


120 ± 20 


Table A3. Spitzer /IRAC infrared photometry used for SED fitting. All fluxes are given in pjy and are corrected for lensing. Upper 
limits are given at the 3cr level. 


Name 

IRAC 

IRAC 

IRAC 

IRAC 


3.6/rm 

4.5/xm 

5.8/xm 

8/im 

MACS0744-system3 

1900 ± 300 

1400 ± 300 



MACS 1149-arcAl. 1 

2000 ± 400 

1500 ± 300 



MACS0451-system7 

360 ± 70 

380 ± 70 



A1413-arc2.1a 

480 ± 90 

490 ± 90 



A1689-arc2.1 

< 340 

< 270 

< 700 

< 700 

A1689-arcl.2 

140 ± 30 

140 ± 30 

< 700 

< 640 





